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Experimental data on the duration of individual fades in a 4-GHz 
radio signal are used to derive the rates of change in amplitude and 
phase of a signal propagated through a 30-mile multipath medium. A 
by-product is a relation between the rate of change in phase and the 
effective vertical velocity of the atmospheric irregularities that can be 
considered to be the cause of the multipath phenomena. The derived 
distribution of velocities is then used to predict the variations in fade 
duration over a wide range of frequency and path length. 

i. introdtjction 

Multipath fading on line-of-sight paths causes substantial variations 
in the amplitude and phase of the received signal. Theoretically, if the 
response to an impulse function were known with sufficient accuracy, 
all magnitudes and rates of change could be computed. In practice, 
this possibility is somewhat of an illusion because many different 
impulse responses occur on any given path, with significant diurnal 
and seasonal variations. Moreover, the desired accuracy is generally 
not achieved, because either the required impulse function needs to 
bo only one cycle of the radio frequency or extensive correlation 
techniques are needed to achieve a comparable result. 

The quantitative estimates given below are based on elementary 
theory and on the measured distribution of fade durations. The results 
agree with (or at least do not contradict) the available experimental 
data. 

II. GENERAL CONSIDERATIONS 

The signal received through a multipath medium can be represented 
by the sum of the principal wave plus the delayed components (a m e ,A '"): 
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E= A3 ,tm +' mU (l + Eay (w+Au,A '") (1) 

= Avf im * MI 

where a,„ and Af,„ represent the magnitude and relative time delay of 
the with "echo." 

The effect of frequency deviation, Am/2tt, is discussed in a later 
section on the variation in net loss within the modulation bandwidth. 
For the single-frequency case (A«> = 0), it is convenient to define 
R, L, (j>, and y as follows: 

v = 1 + ffe' , * +( "- ,) " = Le- iy , (2) 

R sin 



tan 7 = 



1 — R cos <£ 



where R and </>, respectively, are the instantaneous amplitude and 
residual phase of the "composite echo" given by the summation term 
in (1) with Aw = 0. For the deep fades L « 1, | 1 - B | ^ L, | | ^ L, 
and n is an even integer. 

It is shown in Appendix A that 



dL d4> T • / T \ dR 1 

Ldy d4>[, n , ■ rffll 



(3) 

(4) 



and 



where 



= tan (a — /3) 



Ldy 



, , 8ffi7 1 . -i f S1I10 

a = tan I 7 1 = tan I 7 5 

\coa 7 — L/ Vcos <t> — R 

_ 1 / tan 7 \ 

= tan ( zrf I & 7, 

\1 — L/cos 7/ 

At the bottom of a fade dL/dt = and L dy/dt has a maximum 
value of 

d<j) R 

dt (cos 7 — L) ' 
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conversely, dL/dt has a maximum value of 

37 when dy = 0. 

at sm 7 

Frequency selective fading requires that d<j>/dt » dR/dt(fi « 1), 
but dR/dt may not be entirely negligible at the very deep fade levels. 
This observation is consistent with physical models of "reflecting" 
layers that are slowly rising or falling with time. The physical model 
suggests that during any particular fade d<f>/dt and dR/dt can be 
constant, or at least they vary much more slowly with time than the 
"non-linear" parameters L, y, dL/dt, and dy/dt. This information is 
useful in examining the detailed wave interference characteristics 
during a few seconds or minutes in time. On the other hand, the 
process is nonstationary and the average fading characteristics taken 
over many fades depend on the long term average values of d<f>/dt and 
dR/dt, and hence on meteorological conditions. 

Since the variations in signal level are caused by variable metero- 
logical conditions along the path, experimental data are needed to 
obtain quantitative results. In principle, it should be possible to pre- 
dict the radio results from meteorological theory, if the meterological 
parameters were known with sufficient accuracy. In practice, this 
indirect approach is not feasible and direct radio measurements are 
required. Data on the number and durations of deep fades will be used 
to estimate the probability distribution of cos y and d^/dt for use in 
(3) and (4). 

III. NUMBER OF FADES VS DEPTH OF FADE 

Experimental data show that 20-dB fades occur about ten times 
more often than 40-dB fades and on the average last about ten times 
longer. In other words, both the number of fades and their average 
duration are proportional to L. 1 Of all the fades reaching a given level 
L, , 30 percent will "bottom out" at or before 0.7L, , 50 percent will not 
go deeper than 0.5L, , and only 10 percent will reach or go deeper than 
O.lLi . Thus the probability can be written as 

Prob (^ § X) = X 

where X varies uniformly from ^ X ^ 1. 
It can be seen from (2) that 

L = (1 - R cos</>)\/l+tan 2 7 = 2 = R cos * = Rsin * 

cos 7 sin 7 
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and <f> ^ L. At the bottom of a fade, </> -» 0, dL = 0, 7 = 0, and 

L min « 1 - R . 

It follows that at any level L x < 0.1 the ratio 

L mia I 1 - Ro \ 
L x \1 — iti cos0i/ 

^ cos 71 for -Ri tt R . (5) 

Thus the experimentally determined cumulative distribution for 
Anin/^i also applies to cos y, when di?/di « d<j>/dt; that is 

Prob (cos y ^ X) = X 

and the average values are 



^^ cost = f XdX = |. 

i/l •'o 



The corresponding average value of sin 7, needed in a later section, 
is given by 



[Vi-xUx-l 



IV. DURATION OP FADES 



The average duration of a 40-dB fade (L = 0.01) at 4 GHz on a 
30-mile path is about 3 or 4 seconds, and the average duration of a 
20-dB fade is about 10 times longer. The results of an extensive pro- 
pagation experiment have shown that the distribution of fade dura- 
tion t follows a log normal distribution with a standard deviation of 
about 10 log(T 1(i /T 5 o) = 5.6 or In (t 16 /t 50 ) = 1.3. Quantitatively, these 
results are summarized as follows for 4 GHz on a 30-mile line-of-sight 
path : 

99 percent of the fades are shorter than 4000 L seconds, 

70 percent of the fades are shorter than 400 L (average duration), 

50 percent of the fades are shorter than 200 L (median duration). 

These data are plotted on Fig. 1 and it is assumed that fade durations 
shorter than the median are given by the extension of this log normal 
distribution with the same standard deviation.* 



* Recent unpublished data on the duration of rapid fades indicate that the 
standard deviation may be somewhat larger than the "5.6 dB" used in this paper, 
which means that the lines on Fig. 1 might be more nearly vertical than shown. 
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Fig. 1 — Duration of fading at 4 GHz on a 30-mile path. 



The voltage-time diagram .shown on Fig. 1 defines L, L min , <f> , and 
the fade duration t for a fixed frequency. From the diagram it can be 
seen that for dR/dt = the fade duration t is given by 



, 2 /•■'■• 



and since the average value 



y,/R) 



64 



2 Li sin 7, 
d$ 
eft 



B 



(6a) 



sin 7, = 



4 ' 



d(f> _ irL\ 
dt == 2tR 

Substituting into (3) and (4), the average values arc 



dL_ = d(h\L)_ = j*L 
Ldt ~ dt StR' 



(6b) 



(7a) 
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dry _tt_ (1_^_2^) ^ _JL 
dt ~~ 2tR 2 ~4t# 

It follows that, since 20 log L = 8.7 In L dB, 



for L « 1 . 



(7b) 



d7 2 d(ln L) = 2dB/s = dB/s , g) 

dt~ir dt 8.7t " 13.6' 

Figure 2 shows the distribution of rate of change of phase dy/dt and 
the rate of change of L in dB/second, based on (7a) and (7b) and 

on the distribution of fade durations shown on Fig. 1. 

It will be noted that if y were uniformly distributed (sin y = cos 7) 
it would be expected that 

d^ 



dB/s = 8.7 



dt 



but this assumption is not consistent with the experimental result that 
the number of fades is proportional to L. 

In a similar manner, when d<f>/dt = (|^| is a constant < Z/ min ) 
the fade duration might be represented by 



-aC*-a«-« 



dt 



dt 



or 
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Fig. 2— Rates of change in amplitude and phase during fading at 4 GHz on a 
30-mile path. 
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dR 2L 



la"\ 



R t - R„ , (0 . 



dt t _ L cos <f> 

There is no obvious way to relate a change in R to a change in mete- 
orological conditions and there is doubt that it is important to do so. 
The frequency selective nature of fading requires that dxf>/dt be con- 
trolling most of the time. In addition, essentially all fast "flat" fading 
can be explained as a phase interference effect with a delay of one 
or three half-cycles as well as by a supposed change in R. Moreover, 
the experimental results that the number of fades and their average 
durations are both proportional to L are characteristic of a Rayleigh 
type distribution, which is usually explained on the basis of near- 
random phase. Consequently, it is assumed in the remainder of the 
paper that 

di dR 
dt y> dt ' 

V. VERTICAL VELOCITY OF ATMOSPHERIC IRREGULARITIES 

Although dL/dt and dy/dt are important parameters in the design 
of communication systems, the search for a better physical under- 
standing of the fading mechanism needs to begin with the less com- 
plicated parameters R, <f>, and d<f>/dt of the composite echo given by the 
bracketed term in (1). It may be seen from (4) that for dR/dt = 0, 

d<t> L dy 1 



dt dt (cos 7 — L) 

which on the average is 



dt dt (lu; 

In principle, the rate of change in phase d<j>/dt might also be 

obtained by either: (i) direct measurement of phase, or (ii) compu- 
tation from known or assumed values for vertical motion in the 
atmosphere. For example, it is reasonable to expect that 

S = — s2 an 

dt H n - #„_, W 

or conversely, 

. ^ [U - I). \;i, 



IT J dt 
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where df/dt is the rate of change in phase caused by an echo source 
moving vertically through a thickness (H n — #„-i) at a velocity v. H„ 
is the height of the nth Fresnel zone and H„ - H„. x is the thickness of 
the nth Fresnel zone. The relation between d<f>/dt and the effective 
velocity v is shown on Fig. 3 and the probability of exceeding specified 
values of d<f>/dt has been obtained from (10) and from Fig. 2. It ap- 
pears that a vertical motion of only a few centimeters per second (a 
few feet per minute) is sufficient to explain the fade duration measure- 
ments and the values of dy/dt and dB/s derived therefrom. An al- 
ternative derivation of (11) is given in Appendix B. 

VI. FADE DURATION AT OTHER FREQUENCIES AND PATH LENGTHS 

The distribution of velocities obtained in the preceding paragraph is 
expected to be a meteorological phenomena independent of radio 
frequency. It should be possible to use this information to predict the 
fade duration at other frequencies and path lengths. 

Substituting (6b) into (11), and using equation (17), it follows 
that the 

Average Fade Duration = - — (H n — #»-0 = £j -\j~^ ( 12 ) 
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Fig. 3 — Change in phase at 4 GHz vs vertical velocity on a 30-mile path. 
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where D is the path length and X is the wavelength. Equation (12) 
indicates that the fade durations vary as the \A D, if the effective 
value of n is essentially independent of frequency for the average fade. 
The preceding information and assumptions lead to the following 
estimates for the average fade duration at various frequencies, when 
the path length D is measured in km. 



Frequency 


Average Fade Duration 
in Seconds 


160 MHz 


2000 L V| 


4 GHz 


400 L \! 


6 GHz 


320 L V| 


6 X 10 5 GHz 

(visible light) 


*4 



The experimental results at 4 and 6 GHz are not accurate enough to 
establish the A 1/2 variation over this narrow range, but at lower fre- 
quencies the fades, when they occur, last longer than at 4 and 6 GHz 
so the fading rate is not independent of A. 

An extrapolation of the radio data to optical frequencies indicates 
that the average fading rate on a 5-km path through the atmosphere 
would be of the order of one or two cycles per second with much more 
rapid fluctuations for a small percentage of the time. In principle, 
scintillations are rapid multipath phenomena. Although differences, 
such as the effect of water vapor, antenna beamwidth, etc., need to 
be included, the fading rate must vary with the wavelength as A p , 
where p is a positive fraction less than one. A variation of as much 
as A 1 is unrealistic because it would predict that the twinkling of the 
stars would ordinarily be above the flicker rate acceptable by the 
human eye. A consideration of both radio and optical phenomena 
requires a variation in fade duration with frequency close to A 1/2 (at 
least in the range between A 1/:i and A 2/3 ), so the elementary result 
shown in (12) is a reasonable working hypothesis. At optical fre- 
quencies the Fresnel zones are more nearly comparable in size to the 
atmospheric irregularities and both horizontal and vertical motion of 
atmospheric irregularities produce significant scintillations. 



2048 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

VII. VARIATIONS IN NET LOSS WITH FREQUENCY 

The principal interest thus far has been in the characteristics of 
fading at a single fixed frequency. Even if an automatic equalizer 
could be built to compensate perfectly for the effects of fading at any 
reference frequency, the resulting equalization would be less than per- 
fect at frequencies other than the reference frequency. 2 " 4 This type of 
distortion results in variations in net loss (deviations from a constant 
output level) over the modulation bandwidth. By the use of the con- 
cept of Fresnel zones and some simplifying assumptions, it is shown in 
Appendix C that the variations in net loss relative to a fixed level 
at the carrier or pilot frequency is given approximately by 



20 log 






(13) 



M>1 (! + (¥] )• "" 



where 



o„n — , Ttn c ,i . 

20 log L is the depth of fade in decibels (L « 1) and That is the effec- 
tive number of half-wavelengths in the relative delay of the strongest 
delayed component (or composite "echo") that is causing the fade 
at the reference frequency. The corresponding variations in net loss 
are shown in Table I and on Fig. 4. Table I indicates that most of 
the deep fades are caused by "echo" delays of 1, 3, or 5 half-wave- 
lengths (rieff ^ 5) , because experience has shown that good diversity 
requires a frequency separation greater than 1 percent and that accept- 
able transmission quality can almost always be obtained with band- 
widths of 0.1 percent or less. The corresponding values of A/// for a 
20-dB fade are ten times larger than shown in Table I. 

A comparison between theory and experiment is shown on Fig. 4 
for the case of 40-dB fades. This figure also suggests that the median 
40-dB fade corresponds to 1 ^ n e!( ^ 5, but that 10 percent of the 
40-dB fades correspond to delays of many half-wavelengths. The curve 
for n = 100 represents the approximate upper boundary set by 40-dB 
antenna patterns at 4 GHz on a 30-mile path. The corresponding 
theoretical curves for a 20-dB fade would be shifted one decade to the 
right. The crossover of the experimental 40-dB fade data and the 
n = 1 curve shown on the figure is not possible with a single echo, 
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Table I — Variations in Net Lo.ss 





Ne( Loss 

Variation 

(dB) 


Quality 


^for 40-dB fade 


«eff 
L 


Wofl = 1 
(%) 


«eff = 3 

(%) 


tied = 5 

(%) 


0.1 
0.:5 

1 .0 

3.0 

10.0 


0.04 
0.4 
3.0 
10.0 

20.0 


Excellent transmission 
Acceptable transmission 
Poor transmission 
Fair diversity 
Good diversity 


0.03 

0.1 
0.3 
1.0 
3.0 


0.01 

0.03 

0.1 

0.3 

1.0 


0.006 

0.02 

0.06 

0.2 

0.6 



and is an indication that more than two components are present much 
of the time. The "experimental 40-dB fade" frequently represents 
the combination of two or more echoes whose effects add to 40 dB; 
in this case, the limiting theoretical curve for n = 1 is shifted to the 
right because the principal component of the composite echo, taken 
by itself, represents a fade of less than 40 dB. 

The agreement between theory and experiment is considered satis- 
factory and as evidence that most of the fades are caused by strong 
components delayed by only a few half-wavelengths. It has been 
suggested that a better fit with the shape of the experimental data 
exceeded 10 percent and 50 percent of the time can be obtained if the 
i = v— 1 in (13) could be omitted; such an assumption is rejected 




Fig. 4— Variations in net loss during 40-dB fades at 4 GHz on a 30-mile path. 
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because it seems contrary to logic and because it would create dif- 
ficulties in explaining the data exceeded by, say, 90 percent of the fades. 

VIII. CONCLUSION 

Experimental data on the number and duration of fading have been 
used to estimate the rate of change in phase and the rate of change in 
amplitude of an electromagnetic wave received through a multipath 
medium. The elementary theory used herein achieves quantitative 
results, primarily because it explicitly includes the first-order effects of 
phase associated with the use of Fresnel zones. Alternative methods 
starting from the more fundamental wave equations are more sophis- 
ticated and may be more exact in a mathematical sense, but their 
inherent generality does not ordinarily lead to results that can be 
compared with experiment, unless somewhat arbitrary "constants" or 
correlation functions are assumed. 

APPENDIX A 

The relative magnitude and phase of a multipath signal can, by 
definition, be represented as 

v = 1 +tfe 1 ' , * +(n - I) " = Le~ iy 
where n is an even integer and R and <f> are, respectively, the instan- 
taneous magnitude and residual phase of the sum of all the delayed 
components. It follows that 

L = -v/(T - R cos <f>) 2 + (R sin <f>f 
= (1 - R cos <f>) Vl + tan^ 
= V(l - R) 2 + 2#(1 - cos0) 



1 — R cos <t> _ R sin 
cos t sin 7 



By differentiation, 



R sin <f> d<t> — (cos <j> - R) dR 
dL = jr- 

(R cos<£ - R 2 + l-2g cos (ft - 1 + 2R cos4>) dR 
= sin 7 d<}> - — I R 

(\ - R cos 4> - L 2 \ dR 
= sin 7 <** - I Z JR 
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7 , / r \ dR 

= sin y d<f> — (cos y — L) -7Z-- 

In a similar manner, 



dR 
L dy = (cos 7 — L) d<f> + sin y -=r-- 



APPENDIX B 



As shown in several textbooks, the height of the nth Fresnel zone at 
the middle of path length D is given by 

*. - fc - ^' < 15 > 

The corresponding phase delay is 0„ = nv, by definition of Fresnel 
zones. 

In order to restrict the parameter n to integral values, we can gen- 
eralize as follows: 



dd 64 

= n-rr — <f>; — = — -p- 
dt dt 



and 



It follows that 



where v is the effective velocity in a plane perpendicular to the direc- 
tion of propagation. 
From (15) we have 

h„ - #„_, = ivr^ (l - yj^ 1 ) 

^4\n? for w>3 - < 17 > 

By substituting (17) into (16), we have 

rf0 TV 



dt H n -//„_, ' 
which is the same as equation (11) . 
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APPENDIX C 

As noted in (1) the signal received through a multipath medium 
can be represented by the sum of the principal wave plus the delayed 
components (a m e ,A '"'): 

E = Ae «»+»-»(i + £ a m e i{ " + *" )At A (18) 

\ m = l ' 

where a,„ and At m represent the magnitude and relative time delay of 
the mth "echo." The term in brackets can also be resolved into Fresnel 
zones so that 

P-1+ £a n e''< u+au,(n * /u) (19) 

ii- 1 

where n is an integer and a„ is the magnitude of the signal in the 
(n + l)th Fresnel zone. The upper limit N is the maximum number of 
Fresnel zones enclosed by the antenna beamwidth; components with 
longer delays, if present, are relatively unimportant because their 
magnitudes are reduced by the off-beam antenna patterns. The ex- 
pression for v can be separated into the following three terms: 

V = 1 + 2 a n e'" r + Z a n e""(e i& - 1) (20) 

where 8 = (A///)rnr. The first two terms in (20) show the fading 
characteristic at frequency /, which in principle can be corrected by 
a perfect automatic equalizer based on a carrier frequency or pilot 
tone. The third term on the right in (20) shows the relative variation 
(distortion) at frequency A/. 
By definition, let the first two terms be 

1 + X>„e'" T = Lf iy . (21) 

Substituting into (20) and introducing an effective value of 8, we have 

= Le- + U ell [Le- -l-± aj" + ± a^'(f^)] 
= Le- - ."I.,, [l - Le- - ± *^~^ ~ l)] ■ M 
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By this series approximation and a suitable choice for the effective 
value of 8 (S ett ), the last term on the right can be made negligible 
during fading conditions (I <5C 1),* so that 



vttLe- iy [l - 



i8,.,,e' 



(23) 



A "perfect" equalizer (at any given fixed frequency) would leave a 
residual distortion of 



-1-* 



L 



(24) 



The phase angle y is zero at the bottom of the fade, is ±45 degrees 
at 3 dB above the bottom, and varies over the range - tt/2 < y < tt/2. 
In the region of good transmission (8 ett /L « 1) it is sufficiently ac- 
curate to approximate the net loss variations as 



20 log v' = 20 log 



1 - 



L 



= 10 log 



1 + 



L 



(25) 



While both S eff and L arc individually less than unity (and usually 
<0.1) their ratio can be greater than unity, and in that case the net 
loss variations calculated from (25) need to be interpreted as the 
median value of a distribution of net loss variations. 
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* An associated assumption is that the maximum value of 5 is restricted to 



5„,ax = ^A'tT < 1. 



